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A formulation of the Ginzburg-Landau- Wilson version of the partition function of a system with 
a continuously varying order parameter as a transfer matrix calculation allows for the application 
of methods based on the Density Matrix Renormalization Group (DMRG) to the calculation of the 
free energy of the 0(1) model. The essence of both the mapping and the DMRG calculation is laid 
out, along with results that validate this strategy. This method forms the basis of a unified approach 
to the crossover from three to two dimensions in an O(l) system with a slab-like geometry. 



The formulation utilized to analyze the critical point of the O(l) system in three and more dimensions differs 
in fundamental respects from the formulation that has proven most fruitful in two and fewer dimensions. In the 
former case, highly accurate results |Q,|J have been obtained with the use of the Ginzburg-Landau- Wilson effective 
Hamiltonian, which is based on a continuously varying order parameter. On the other hand, lower dimensional systems 
are generally modeled with the use of a fixed-length order parameter ||f|]. 

The above state of affairs is acceptable, if not entirely satisfactory, as long as one is not interested in a system 
whose critical point properties possesses both three- and two-dimensional features. Such is the case in the instance 
of dimensional crossover in a system, with a slab-like geometry. Scaling relations and the Renormalization Group 
predict the asymptotic behavior of such a system when the correlation bulk three-dimensional correlation length is 
asymptotically large or small compared to the finite slab width. In addition, finite size scaling indicates the existence 
of a general form for the crossover between those two limits (5|(|. However, the calculation of the specific form of the 
crossover function requires a unified theoretical approach to the critical point properties of a slab-like system in both 
the three- and two-dimensional regimes. 

The calculation of the crossover function, and of the thermodynamic quantities that can be derived from it, has 
been carried out in a few cases. For example, crossover from three to two dimensions has been worked out in the 
mean spherical model HQ and in its mathematical equivalent, Bose-Einstein condensation. Crossover from d to d— 1 
dimensions can also be determined in an 0(1), or Ising-like, system when d > 5 ||]. Here, both systems possess 
asymptotically mean-field critical behavior. 

In addition to the above exact determinations of dimensional crossover, a field-theoretically motivated approach has 
been formulated based on the notion of an "environmentally friendly" renormalization group PJlOfl. This approach 
yields explicit (and relatively simple) crossover functions for thermodynamic quantities. While the functions are nec- 
essarily approximate, they are consistent with the key expectations that arise from finite-size scaling. Unfortunately, 
existing literature provides no simple way to improve the crossover functions that arise from the environmentally 
friendly renormalization group, so as to produce predictions that can be profitably compared with experimental data. 

We have developed an approach that holds the promise of producing just such an improvement. This approach is 
based on a reformulation of the partition function of a lattice based system with a continuous order parameter in terms 
of transfer matrices which allows for the calculation of the partition function either with the use of field-theoretic 
techniques or with the use of calculational devices that have been formulated to speed the numerical evaluation of 
quantum mechanical density matrices in one dimension, a problem that is mathematically equivalent to the calculation 
of transfer matrices in two dimensions. 

We start with the Ginzburg-Landau- Wilson effective Hamiltonian of an O(l) system on a lattice in three dimensions. 
If Si is the value of the continuous spin variable at the i th lattice site, then this effective Hamiltonian has the form 



The first sum in Eq. (gj) is over nearest neighbor pairs. High order perturbation theory coupled with the Renormal- 
ization Group has been shown to yield very accurate results for the three-dimensional critical behavior of the system 
described by the effective Hamiltonian ([!]). 

In two dimensions, the Ising model, which is based on fixed length spins, has been successfully investigated with 
the use of a variety of techniques. An important subset of the approaches is based on an analysis of the transfer 
matrix coupling the state in 2 Ar -dimensional space of a column of spins to the state of a neighboring spin column 
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p|. Onsager's exact solution [Q for the Ising model partition function consisted of the extraction of the largest 
eigenvalue of this transfer matrix. In the absence of techniques leading to exact solutions, one has recourse to other 
methods for the determination of the largest eigenvalue of a two-dimensional transfer matrix. In particular, there is a 
set of approaches developed to perform highly accurate numerical evaluation of the eigenvalues of the density matrix 
of a one-dimensional quantum mechanical system. The generic term for this collection of calculational devices is the 
Density Matrix Renormalization Group (DMRG) , ]l5| . There are cases in which the extremely high accuracy of 
this method has led to important results. 

In the context of statistical mechanics, the DMRG has an analog in the form of the Transfer Matrix Renormalization 
Group (TMRG) [[l3| . The TMRG has been shown to lead to excellent results for the thermodynamic functions of the 
two dimensional Ising model []l3| , the Q-state Potts model [[l6) , |l7j and confinement effects in the presence of gravity 
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We have been able to demonstrate that this method can be applied to the evaluation of the partition function for the 
continuous-spin effective Hamiltonian (jl]) . We start with the one-dimensional version of the Ginzburg-Landau- Wilson 
partition function 



[ J dsi exp 



Si+l, 



2 S * + 4* 



(2) 



This partition function can also be written as the inner product of a sequence of transfer matrices T(s,, Sj+x), where 
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T(si, sj) = exp 

Eigenfunctions, ipk(s), of this transfer matrix are defined by the equation 

ds'T(s,s')ip k (s')=\ k ip k (s) (4) 

We will list the eigenvalues in order of magnitude. The largest eigenvalue will be Ao, the next largest Ai and so on. 
Given that the transfer matrix is real and symmetric, the eigenfunctions can also be written as real, normalized and 
orthogonal, in that 
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ipk{s)ipi(s)ds = 5 k ,i 
The one-dimensional transfer matrix is then written as 



(6) 



Given the orthnormality of the eigenfunctions, the one-dimensional partition function has the form 



(7) 



and is, thus dominated by the contribution associated with the largest eigenvalue, Ao, of the transfer matrix. 

Now, the d-dimensional Ising model can be written in terms of the eigenfunctions and eigenvalues of the one- 
dimensional Ising model by associating neighboring pairs of spin variables in terms of a "bond" transfer matrix, 
Tb(s{,Sj), having the form 
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The partition function is then a sum over bond indices, k, while vertices contain integrations over the spin variable s 
of combinations of transfer matrix eigenfunctions: 
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While the eigenfunctions, ipk(s) are orthonormal, there is no reason to believe that this will lead to significant 
restrictions over the allowed indices in the integration of a product of 2d of the ?/>fc's when d > 2. One important 
simplification arises from the fact that a tpk(s) will be either even or odd in the spin variable s. In fact, the largest 
eigenvalue is associated with an even-parity eigenfunction and the next largest with an odd-parity eigenfunction. 

The partition function of the two-dimensional continuous-spin variable version of the O(l) model was calculated 
with the use of the formulation above in terms of transfer matrix eigenfunctions, the basis set being truncated at the 
eigenfunctions having the largest few eigenvalues. The fourth order coupling constant, u, in (|l|) was chosen to be 
equal to 0.23, guided by a best fit between the Ginzburg-Landau Wilon effective Hamiltonian and order parameter 
distributions in three dimensions as obtained by simulation |2(| . These eigenfunctions were used to construct the basis 
set of the two-dimensional transfer matrix connecting one row of continuous spins to the next. This two-dimensional 
transfer matrix was then reduced with the use of an adaptation of the TMRG. The free energy is extracted from the 
largest eigenvalue of the transfer matrix via the connection T = —ksT In Ao. Because we are interested in the behavior 
of the partition function in the immediate vicinity of the critical point, temperature-dependent factors multiplying 
the log of the largest eigenvalue are set equal to a constant. 

This version of the free energy leads immediately to the specific heat through differentiation with respect to the 
"bare" reduced temperature, r. It is also possible to introduce a symmetry-breaking field, h, appearing in the 
effective Hamiltonian through the term — hsi. The remanent magnetization and the isothermal susceptibility can 
be calculated from the free energy by taking, respectively, first and second derivatives with respect to the symmetry- 
breaking field. 

The spontaneous magnetization of the system can be calculated from the spin operator S2 



where ^0 is the eigenfunction of the transfer matrix with the largest eigenvalue and the spin operator S2 is the Ising 
spin operator. We plot the value of ([h]) on a log-log grid to illuminate the power-law behavior (see Figure |l|). The 
value of the critical exponent (3 is calculated to be 0.1229 ± 0.0006 which is consistent with the exact value for the 
d = 2 Ising model of 1/8 Q. 

Differentiating the free energy twice with respect to the reduced temperature r, we obtain the specific heat. We 
display our results for this quantity as a function of r in Figure Q Further analysis of these data reveal a logarithmic 
dependency both above and below the critical temperature, in other words a = at = 0, as expected from the Ising 
model. 

We have also examined the ratio of specific heat amplitudes A + /A- where we fit to the function 



For an O(l) system in two dimensions, this value is expected to be unity. The value calculated from the data is 0.94. 

As mentioned in the previous section, by adding a symmetry-breaking term — hsi, we are able to take derivatives 
with respect to h to obtain the spontaneous magnetization (independently from the (^ol^l^o) method) and the 
magnetic susceptibility xt- Determination of the magnetization using this approach gives (3 — 0.1230 ± 0.0009. 

The calculation of \t requires special consideration as it is dependent on the size of the "grid" Ah used for 
the numerical differentiation. Best results for the critical exponent above and below the critical temperature are: 
7 = 1.74 ± 0.01 and 7/ = 1.33 ± 0.03. The exponents do appear to converge to the expected value of 7/4 as the grid 
size dh is reduced (see Figure ||) . 

Finally we comment on the validity of the assumption that the retention of a small number of states in the 
eigenfunction expansion of the transfer matrix (^|) suffices to ensure an accurate calculation of critical point properties. 
The previous results were obtained by using only two states in the expansion (|^), i.e. n s = 2. To empirically investigate 
the inclusion of more states, the calculation of the specific heat was repeated with n s = 4. The results are effectively 
the same within error bounds - a = ai = and the amplitude ratio is 0.95. 

There is another perspective from which to explore the effect of truncating the Ginzburg-Landau basis. The lowest 
order effect of including higher states can be realized by allowing a single bond in the lattice to be represented by the 
full sum of states in the expansion. We will refer to this as an "enhanced" bond. All other bonds are expanded on 
only the first two (n s = 2) eigenstates. We can then observe the effect of this enhancement on the free energy density 
as we increase n s . Examining Figure H it is clear that the effect of enhancement is negligible above n s = 4. 
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FIG. 1. The spontaneous magnetization (calculated as the expectation value in the state |^o)) as a function of (r c — r). 
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FIG. 2. The specific heat as a function of r. 
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FIG. 3. The susceptibility exponent 7 above and below the critical temperature as a function of the grid spacing dh. The 
horizontal line denotes the 2D Ising exponent of 7 = 1.75. 
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FIG. 4. The free energy of a chain with a single enhanced bond as a function of n 3 (i.e. the number of states kept in the 
expansion of the Ginzburg-Landau transfer matrix) . The inset shows the detail of n s — 4 and higher. 
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